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ABSTRACT 

The  SDIO  Hi9b  Altitude  Balloon  Experinent  (HABE) 
required  a  rapid  controlled  descent  of  a  large, 
heavy  payload  iron  a  balloon  platform.  An 
existing  viscous  brake  concept  was  chosen  for 
the  job,  but  the  short  tine  and  space  available 
coupled  to  the  heavy  payload  resulted  in  very 
high  power  densities.  Existing  rules  of  thunb 
for  viscous  brakes  were  insufficient  for  the 
design.  Viscous  brake  pbenonenology  and  the 
particular  physical  constraints  were  coded  into 
a  nodern  continuous  simulation  applications 
package.  After  investigating  both  full  dynamic 
and  quasi-steady  models,  it  was  decided  that  the 
latter  was  adequate  for  all  but  the  first  and 
last  small  fraction  of  the  descent.  The  model 
includes  temperature  dependent  viscosity  and  a 
beat  transfer  model.  The  model  been  used 
extensively  to  investigate  the  design  space  and 
to  develop  test  plans. 

BACKGROUND 

A  unique  ascent  procedure  has  been  proposed  for 
the  HABE  experiment.  At  liftoff,  all  of  the 
helium  will  be  contained  in  a  *tow  balloon*.  At 
a  designated  altitude,  the  payload  will  be 
dropped  approximately  500  feet  below  the  tow 
balloon,  allowing  the  deployment  of  the  main 
balloon  and  parachute  which  are  attached  to  the 
tow  and  stored  at  the  top  of  the  payload 
section.  See  Figure  1  for  a  representation  of 
the  reel  down  sequence. 

The  design  of  the  system  is  constrained  by  the 
following:  full  reel  down  must  by  accomplished 
within  3.5  minutes  to  insure  that  the  marn 
balloon  is  fully  deployed  before  gas  begins  to 
flow  from  the  tow  balloon  to  the  main  balloon; 
at  the  end  of  the  SOOft  reel  down,  separation 
velocity  must  be  low  enough  to  avoid  excessive 
stress  on  the  main  balloon  and  parachute  which 
must  provide  the  final  breaking  force;  the  fluid 
temperature  must  not  exceed  450*r  to  avoid 
thermal  breakdown.  A  desired  capability  is  that 
the  system  should  stay  within  these  constraints 
for  initial  conditions  within  the  range  from  0*F 
to  lOOT. 

The  Aerospace  Engineering  Division  of  the 
Phillips  Laboratory  has  had  extensive  experience 
using  reel  down  systems  with  integrated  viscous 
brakes  (Wagner  and  Doherty  1958).  The 
requirements  for  the  H>.BE  mission  are  much  more 
severe  in  payload  weight  and  time  available  for 
the  deployment.  This  causes  very  high  brake 
power  loads,  resulting  in  significant  beating  of 
the  brake  fluid  which,  in  turn,  substantially 
lowers  fluid  viscosity.  If  the  brake  fluid 
viscosity  is  too  low,  thermal  runaway  could 
result.  If  it  is  too  high,  reel  down  will  take 
too  long.  A  mathematical  simulation  of  the_ 
deployment  was  required  to  assist  in  the  design 
and  test  of  the  reel  down  system. 


XUESBI 

The  motion  of  the  system  is  determined  by  applying 
Newton's  2nd  Law  to  each  point  mass: 


Where:  m  ■  mass 

K  B  position  vector  to  center  of  mass 
r,  «  the  i**  force  vector 

Since  reel  down  would  occur  at  some  time  after 
liftoff,  it  is  assumed  that  the  system  would  be 
traveling  horizontally  at  the  speed  of  the  wind, 
and  that  any  horizontal  relative  motion  between 
balloon  and  payload  would  be  damped  out.  This 
reduces  the  scope  of  the  analysis  to  the  vertical 
coordinate  only. 

A  free  body  diagram  of  each  mass  is  shown  in  Figure 
2.  The  equation  of  motion  for  each  mass  is: 

•  B  ~  P.  -  T  ~  (2) 

at* 


■  r  -  (w,  -  dw)  (3) 

o  t 


Where:  T  «  Tension  in  the  connection 

w  »  Weight 
B  “  Buoyancy 
D  °  Drag 
Subscripts; 

1  •  Balloon  and  Helium 
3  Payload 


courLmoN 
Of  acErDowN 


Figure  1.  The  reeldown  sequence. 
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and  where  A  y  the  weight  of  the  main  balloon  and 
parachute  that  has  been  extracted  from  the 
payload  section.  This  added  weight  on  the  upper 
nass  should  include  the  additional  force 
required  to  accelerate  the  main  balloon  and 
parachute  from  the  velocity  of  the  payload  to 
the  velocity  of  the  upper  mass.  This  tension  is 
equal  to  the  instantaneous  rate  of  chanqe  of 
momentum  of  those  items  which  is  negligible  for 
this  problem. 


Figure  3.  Isometric  of  brake  and  reels. 

As  stated  earlier,  the  HABE  brake  will  be  absorbing 
a  lot  of  power,  which  will  cause  fluid  heating. 

The  temperature  dependence  of  viscosity  for  many 
liquids  is  known  to  follow  the  following 
relationship: 


V  •  A  V,  exp  (B/T) 


(6) 


Figure  2.  Free  Body  Diagram  of  balloon  and 
payload. 

Similarly,  it  is  assumed  that  the  combined 
tension  in  the  cables,  r,  is  equal  at  both  ends 
which  implies  that  the  cable  mass  is  negligible. 
The  separation  or  reel  down  speed  is  controlled 
by  the  tension  in  the  line,  which,  in  turn,  is 
controlled  by  the  viscous  brake.  The  brake  is 
on  the  shaft  connecting  the  two  reels  as  shown 
in  Figure  3.  When  Newton's  Second  Law  is 
applied  to  a  system  rotating  about  a  fixed  axis 
the  result  is: 

where:  J  ~  moment  of  inertia  of  the  rotating 
system  (slug  ft*) 

0  >  angle  of  rotation  (rad) 

K  •  radius  to  cable  on  reel  (ft) 

“  moment  caused  by  brake  (lb  ft) 

The  brake  moment  (Wagner  and  Doherty  1958 )is: 


Where  v.  is  the  nominal  viscosity  for  the  fluid  and 
A  and  B  are  constants  for  the  material.  For  the  CC 
VISCASIL*  fluids  used  in  the  HABE  brakes,  the 
constants  were  found  to  be  0.003  and  1740 
respectively  (when  r  is  measured  in  K). 

The  reel  down  is  accomplished  by  spooling  DuPont 
KEVLAR*  cable  from  two  straight  sided  spools. 
Assuming  a  continuous  reduction  of  radius  with 
length,  the  instantaneous  radius  is; 

R  •  -  (rI  -  Kt)l/L  (7) 

where:  1  is  the  instantaneous  length 
L  is  the  total  length 

and  subscript  o  refers  to  the  outer  extent  of  the 
cable  radius  when  fully  wound  on  the  spool,  and  i 
is  the  radius  when  fully  deployed. 

The  rate  of  change  of  length  of  cable  is  related  to 
the  turning  of  the  cable  by  the  rotational  formula: 


•  fl  n  V  p  o  (ri-ri)/d  (  5  ) 


Where:  r  ■  kinetic  viscosity  (ft’/s) 

p~  brake  fluid  density  (slugs/ft*) 
w  “  angular  velocity,  same  as  dS/dt 
(rad/s) 

a  ■■  nusd>er  of  rotors  in  the  brake 
r  "  outer  or  inner  radius  of  active 
brake  surface  (ft) 

d  ■  distance  between  disk  rotors  and 
stators  (ft) 


If  the  cable  is  assumed  to  be  have  infinite 
stiffness,  then  the  separation  of  the  two  masses  is 
exactly  deteridned  by  the  amount  of  reel  down, 
which  reduces  the  degrees  of  freedom  of  the  problem 
to  2.  This  seems  especially  justified  since  DuPont 
KEVLAR*  is  a  very  stiff  material.  With  infinite 
stiffness,  the  length  rolled  off  the  spool,  I,  is 
directly  related  to  the  distance  between  the  two 
masses  by  the  equation: 
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1  -  Yi  -  y. 


T  R  d 


(15) 


(9) 


u 


/I  n  V  p  (rj  -  ti) 


The  approach  used  to  Bolve  the  differential 
equation  set  was  to  first  solve  Equation  4  for 
the  tension,  r,  and  substitute  the  result  into 
the  two  linear  equations  of  notion.  Equations  2 
and  i.  Next,  Equations  8  and  9  are  used  to 
relate  the  angular  acceleration  to  the  two 
linear  accelerations  as  follows:  Take  the 
second  time  derivative  of  Equation  9  to  obtain: 


d^l  .  tfVi  _  d*y, 
d  t*  d  t’  d  t* 


(10) 


The  time  derivative  of  Equation  8  is: 

cf»i  .  dR  de  ^  -  d»e 
d  t*  d  t* 


(11) 


Combining  Equations  10  and  11,  and  solving  for 
the  angular  acceleration  tern,  we  have: 

.  tfVi  _  dVi  _  dR  ^  f  12) 
d  t’  d  t’  d  t*  dt  '  ' 


Similarly,  the  steady  state  assumption  is  applied 
to  the  equation  of  motion  for  the  bottom  nass 
(Equation  3)  to  yield  the  result  that  the  tension 
is  equal  to  the  weight  of  the  payload  minus  the 
amount  of  balloon  and  payload  that  has  been 
deployed.  Note  that  this  steady  state  solution 
does  not  result  in  a  constant  velocity.  As  the 
cable  pays  out,  T  and  R  decrease,  which  tends  to 
decrease  the  angular  velocity.  The  increasing 
temperature  lowers  viscosity  which  tends  to 
increase  angular  velocity.  The  solution  is 
"steady”  in  that  the  changes  are  slow  enough  that 
the  acceleration  terms  are  not  a  significant  part 
of  the  equations  for  most  of  the  time  period  of 
interest. 

To  complete  the  description  of  the  steady  state 
solution  set,  the  initial  conditions  completely 
determine  the  right  hand  side  of  Equation  IS.  As  4i 
(i.e.,  dO/dt)  is  determined.  Equation  8  is 
integrated  to  determine  the  amount  of  cable  that 
has  been  reeled  out  so  that  the  instantaneous  value 
of  r  and  R  can  be  determined. 

The  instantaneous  brake  power  for  rotating 
equipment  is  the  torque  times  the  angular  velocity: 


Equation  12  is  then  used  to  remove  angular 

acceleration  from  the  linear  motion  equations  to  Ptrua  “To  (16) 

obtain  Equations  13  and  14: 


K 


R’J  dc 


J  dVj 
R*  dt 


B  -  Wi-  AV 


2LR  \Tt)  ~T- 


(13) 
-  D 


m. 


-  d  dVj 
R’J  dt  dt 


-  If,*  AH 


(14) 


d{Rl  -  Ri*)/ce\»  *  Mtru, 
2  L  R  \3t)  R 


The  differential  equation  solution  software, 
ACSL*,  can  not  solve  these  equations  as  written, 
since  both  linear  accelerations  appear  in  each 
equation.  This  is  easily  remedied  by  using 
algebraic  methods  to  obtain  an  equation  for 
either  the  tine  derivative  of  V,  or  V,.  Only 
one  of  them  must  be  isolated  for  ACSL  to  solve 
the  problem. 


This  dissipated  power  appears  as  heat  in  the  brake 
fluid.  As  the  fluid  heats  up,  it  transfers  heat  to 
the  rotors  and  stators  which,  in  turn,  pass  it  on 
to  the  housing  and  the  shaft.  The  heating  period 
is  so  short  that  not  very  much  is  transferred  to 
the  surrounding  air.  The  heatin:)  of  the  brake  was 
investigated  by,  first,  formulating  a  detailed 
thermal  model  of  the  brake,  then  using  the  results 
of  that  model  to  guide  in  the  formulation  of  a 
simpler  model  for  incorporation  in  the  ACSL 
solution.  That  model  was  a  three  node  model  which 
consisted  of  three  equations  of  the  form: 

*0i  ( 17  ) 

Where:  r,  »  temperature  of  the  i*‘  node 

K,i  -  thermal  conductance  from  1'“  to  node 
e,  =  internal  heat  generation  for  the  i“' 
node 

M,  =  mass  of  the  i'‘  node 

c,  “  heat  capacity  of  the  i'*  node 


A  constant  temperature  solution  of  the  equations 
revealed  that  the  system  very  rapidly  converges 
to  the  steady  state  solution  of  the  angular 
velocity  equation.  Equation  4.  The  time  to  get 
to  this  steady  state  solution  is  less  than  one 
second  out  of  a  complete  run  of  almost  200 
seconds,  as  shown  in  Figure  4.  This  is 
fortunate  since  the  steady  state  solution 
requires  much  less  computation,  and  is  much 
easier  to  use  in  studying  the  effects  of 
temperature  and  other  variables. 

Steady  State  Solution: 

The  steady  state  solution  of  Equation  4  is 
obtained  by  setting  the  second  derivative  of  0 
to  zero,  which  yields  the  static  result  that  the 
sum  of  the  moments  is  zero.  That  is,  the  torque 
from  the  reel  down  cables  is  exactly  balanced  by 
the  torque  from  the  brake.  After  solving  for 
angular  velocity,  we  obtain: 


The  nodes  represent  the  following:  1,  the  fluid, 
which  was  the  only  node  to  have  any  internal  heat 
generation,  2,  the  rotors  and  stators,  and  3  the 
housing  and  shaft.  The  values  used  for  thermal 
mass  (M,C,)  were  3,990J/K  for  1,  21,100J/K  for  2, 
and  14,600J/K  for  3  (Hironcr  1992).  The 
conductance  between  1  and  2  was  680K/K,  between  1 
and  3  was  47.5H/K,  and  between  2  and  3  was  9.50W/K. 
The  internal  heat  generation  is  the  power 
dissipation  computed  in  Equation  IS. 

INITIAL  DESIGN  CALCULATIONS 

The  steady  state  program  was  run  with  various 
nominal  fluid  viscosities  to  determine  what  was 
required  for  acceptable  operation  of  the  system. 
With  an  initial  temperature  of  60*F,  and  a  nominal 
viscosity  of  12S,000  centistokes,  the  SOO  foot  line 
deployed  in  200  seconds.  With  both  the  payload 
weight  and  the  radius  of  the  line  on  the  reel 
decreasing  with  time,  the  torque  on  the  brake  is 
monotonically  decreasing  with  time  as  shown  in 
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Figure  5.  If  viscoBir.y  were  constant,  brake 
reel  speed  would  also  decreased,  but,  in  fact, 
the  reel  speed  increases  as  shown  in  Figure  6. 
The  reason  for  this  is  that  the  power 
dissipation  into  the  brake  (Figure  7)  is 
increasing  the  temperature  of  the  brake  fluid 
(Figure  8),  which  is  decreasing  the  viscosity. 
The  problem  is  fully  coupled  in  that  the 
increasing  rollout  velocity  increases  the  power, 
further  increasing  the  heating,  which  further 
reduces  viscosity. 

while  the  above  description  suggests  a  potential 
run-away  situation,  the  final  result  is  fairly 
benign  with  the  proper  brake  fluid  viscosity. 

The  decreasing  torque  finally  balances  the 
effect  of  increasing  temperature  to  slow  down 
the  rotational  speed.  The  decreasing  radius 
also  results  in  decreasing  rollout  velocity  as 
shown  in  Figure  9.  This  particular  case  is 
judged  successful  since  brake  fluid  temperature, 
total  time,  and  final  rollout  velocity  were 
within  the  design  constraints.  This  viscosity 
will  behave  differently  at  different  starting 
temperatures,  as  shown  in  Figure  10. 

PROTOTYPE  TESTING 

At  this  point  it  is  proper  to  point  out  that  the 
model  described  above,  while  sophisticated  for 
standard  brake  design  practice,  contains  many 
assumptions  and  simplifications  which  lead  to 
inaccuracies  in  the  final  results.  Some  of  the 
major  uncertainties  are  the  small  differences  in 
geometry  from  one  rotor-stator  worki'.g  surface 
to  the  next  due  to  manufacturing  tolerances,  the 
actual  viscosity  of  the  fluid  at  the  shear  rates 
of  the  test,  the  accuracy  of  a  single  average 
temperature  in  determining  the  correct  average 
viscosity  of  the  fluid  in  the  brake,  the 
accuracy  of  the  3  node  model  for  predicting  the 
heat  transfer  events  of  the  btaKe,  and  the 
accuracy  of  the  thermal  masses  'and  the  thermal 
conductances  in  the  model.  A  horizontal  brake 
reel  roll  out  test  was  planned  to  insure  that 
the  system  and  components  performed 
satisfactorily  under  loads  expected  for  the 
operational  system  and  to  validate  the  system 
model.  Only  by  this  development  testing  could 
we  have  confidence  of  success  when  the  system  is 
dropped  from  a  balloon. 

Two  trucks  were  used  to  perform  the  horizontal 
brake  reel  test.  One  truck  was  used  to  hold  the 
system  in  place  and  another  truck  pulled  out  the 
Kevlar  Line.  Since  the  trucks  were  available 
for  only  a  limited  amount  of  time,  and  since 
each  test  will  require  at  least  2  hours  for  cool 
down  of  the  brake  and  rewinding  of  the  line  on 
the  spools,  the  probable  number  of  tests  is  from 
2  to  4.  In  order  to  design  the  test  for 
maximum  useful  data,  several  possible  test  runs 
were  simulated  with  the  system  model.  Accurate 
speedometer  and  load  cell  readout  was  available 
in  the  truck  cab. 

The  basic  model  was  modified  to  simulate 
constant  applied  force,  constant  pull  out  speed, 
and  constant  applied  power.  The  constant  force 
resulted  in  a  speed  curve  which  was  similar  to 
the  expected  operational  performance  shown  in 
Figure  9.  The  resulting  power  input  was 
proportional  to  the  speed  curve.  Constant  speed 
operation  resulted  in  high  initial  force  which 
decreased  to  a  shallow  minimum  near  the  midpoint 
of  the  pull,  then  increased  dramatically  near 
the  end  of  the  run.  The  constant  power  runs 
resulted  in  an  increasing,  then  decreasing  speed 
requirement  and  a  decreasing,  then  increasing 
force  input. 

The  best  option  for  calibration  of  the  thermal 


model  was  the  constant  power  run.  It  resulted  in 
the  longest  period  of  nearly  linear  temperature 
rise,  but  this  condition  imposes  the  requirement  of 
either  maintaining  a  variable  speed  or  a  variable 
pulling  force,  which  would  probably  be  difficult 
for  the  driver  to  accomplish.  For  this  reason, 
this  option  is  not  recommended  for  the  test. 

The  constant  speed  tests  had  the  advantage  of  a 
considerable  long  period  of  nearly  constant  load, 
followed  by  a  gradually  increasing  load  toward  the 
end  of  the  run.  This  allowed  a  thorough  mechanical 
test  of  the  system  at  nearly  nominal  loads  followed 
by  a  final  increase  to  proof  test  load.  This 
condition  also  results  in  a  nearly  linear 
temperature  increase  over  a  long  part  of  the 
experiment,  which  would  facilitate  calibration  of 
the  thermal  model.  The  disadvantage  of  this  test 
option  is  that  the  chosen  speed  must  be  low  enough 
to  avoid  damaging  loads  at  the  start  and  end  of  the 
run  which  keeps  average  power  and  total  energy 
levels  lower  than  that  expected  on  the  HABE 
deployment  drops,  which  results  in  small  increases 
in  brake  fluid  temperature. 

The  constant  pull  force  option  has  the  advantage 
that  the  average  power  and  total  energy  levels  can 
be  of  the  same  magnitude  as  the  actual  HABE 
deployments.  This  gets  the  fluid  temperature  up  to 
expected  operational  levels,  thereby  validating  the 
brake  torque  model  throughout  the  operational 
envelope. 

It  is  assumed  that  the  driver  will  have  large 
readout  gages  for  both  force  and  speed  in  the  low 
range  of  interest.  While  the  driver's  ability  to 
maintain  either  a  constant  speed  (with  varj'ing  pull 
force)  or  a  constant  pull  force  (with  varying 
speed)  is  not  assured,  so  long  as  the  actual  speed 
and  force  data  are  known,  the  model  can  be 
calibrated. 

CONCLUSION 

The  continuous  simulation  software  proved 
invaluable  in  both  the  design  and  testing  of  the 
new  brake  assembly.  In  fact,  without  a  simulation 
model  it  might  have  been  impossible  to  arrive  at  a 
satisfactory  design  for  the  brake,  or  even  be  able 
to  interpret  the  test  data. 
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Figure  4.  ComperisoD  of  dyoanic  to  steady  state 
solution  for  a  constant  temperature 


Figure  8.  Brake  fluid  temperature  during  rollout 


